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ABSTRACT 

Using the 8.4m Large Binocular Telescope, we observed six GRB afterglows 
from 2.8 hours to 30.8 days after the burst triggers to systematically probe the late 
time behaviors of afterglows including jet breaks, flares, and supernova bumps. 
We detected five afterglows with Sloan r' magnitudes ranging from 23.0-26.3 
mag. The depth of our observations allows us to extend the temporal baseline 
for measuring jet breaks by another decade in time scale. We detected two jet 
breaks and a third candidate, all of which are not detectable without deep, late 
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time optical observations. In the other three cases, we do not detect the jet 
breaks either because of contamination from the host galaxy light, the presence 
of a supernova bump, or the intrinsic faintness of the optical afterglow. This 
suggests that the basic picture that GRBs arc collimated is still valid and that 
the apparent lack of Swift jet breaks is due to poorly sampled afterglow light 
curves, particularly at late times. 

Subject headings: gamma rays: bursts 

1. Introduction 

It is a basic ingredient in GRB models that GRBs are collimated in jets (e.g., Rhoads 
1999). This assumption affects the estimates of both the total energy of GRBs and the rate 
of GRBs. Depending on the GRB jet model, a beaming correction of ~ 20-500 is needed 
to obtain the intrinsic rate of GRBs (e.g.. Frail et al. 2001; Granot et al. 2002; Zhang et 
al. 2004; Yamazaki et al. 2004; Dai & Zhang 2005), and uncertainties in this correction 
need to be resolved before wc can fully understand GRBs and their progenitors. Before the 
launch of Swift, jet breaks had been observed in many optical light curves (see Zeh et al. 
2006 for a collection of pre- Swift jet breaks) and several bursts had achromatic breaks across 
the optical bands (e.g., GRB 990510, Stanck et al. 1999), consistent with the predictions of 
the jet model. After the launch of Swift (Gehrels et al. 2004), with its rapid localization of 
GRBs and the dedicated on-board XRT instrument, the number of GRBs with jet breaks 
was expected to increase significantly. Instead, the afterglow light curves were found to be 
more complicated, especially in the X-ray band, where multiple breaks or giant X-ray flares 
are present in more than half of the Swift-XBT light curves (Nousek et al. 2006; O'Brien 
et al. 2006). The optical afterglow light curves also show complicated decay patterns both 
in pre- Swift bursts (e.g. Bersier et al. 2003) and Swift bursts (e.g., Stanck et al. 2007). 
With these complications in the afterglow light curves, it can be difficult to identify a jet 
break. As a result, the number of jet breaks detected in the Swift era is rather small (e.g., 
GRB 050525A, Blustin et al. 2006, GRB 060206, Stanek et al. 2007, and GRB 060526, Dai 
et al. 2007, and see Liang et al. 2007 for a comprehensive analysis). Indeed, the absence of 
late time X-ray breaks can be interpreted as a challenge to the basic beaming model (e.g.. 
Burrows & Racusin 2007a) or the fireball model (Dado et al. 2007). 

While the fiares are compfications, the biggest problem for identifying jet breaks is the 
need for well-sampled, long term light curves. While almost all Swift bursts have good 
signal-to-noise ratio (S/N) X-ray light curves at early decay times, most of the jet breaks 
occur at late times when the uncertainties in the XRT light curves increase significantly. 
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In these cases, the claim that the data are consistent with the single power-law decay does 
not distinguish between models, because the error-bars are so large that the data are also 
consistent with the broken power-law of a jet break — jet breaks are "hidden" in low S/N 
light curves or "masked" by additional sources of emission (Curran et al. 2007; Ghirlanda et 
al. 2007; Sato et al. 2007). Shao et al. (2007) have also proposed that jet breaks are hidden 
by X-ray light echos from dust scattering. In the optical bands, the afterglow monitoring 
is generally sparse, in part because of the large number of bursts detected. However, the 
late-time afterglow can be measured accurately using large optical telescopes. In fact, all 
the Swift bursts where jet breaks have been detected have well-sampled optical light curves 
(e.g., Dai et al. 2007). 

Given the need for deep, late-time monitoring of GRBs, we initiated an observing pro- 
gram using the newly built 8.4 meter Large Binocular Telescope (LET, Hill et al. 2006). 
With 20-30 minute exposures, we are able to reach a flux hmit of ~ 25-26 mag in the Sloan 
r' band. This allows us to extend the temporal baselines of optical afterglow monitoring 
by roughly an order of magnitude, and systematically search for the jet breaks and other 
emission components in the optical light curves. 

2. Observations and Data Reduction 

We observed 6 GRBs, GRB 070125 (Hurley et al. 2007), GRB 070311 (Mereghetti et al. 
2007), GRB 070411 (Moretti et al. 2007), GRB070412 (Romano et al. 2007), GRB070419A 
(Stamatikos et al. 2007), GRB070518 (Guidorzi et al. 2007), with the 8.4 m LBT (Hill 
et al. 2006) using the Large Binocular Camera (LBC) Blue CCD camera (Ragazzoni et al. 
2006) during the LBT/LBC Science Demonstration Time period in January-November 2007. 
A journal of the LBT observations is presented in Table 1. We took 5-25 dithered, 200 s 
exposures for each epoch using the SDSS r-band filter. All the images from a given epoch 
were registered to a common astrometric grid and co-added. We ran DAOPHOT/ALLSTAR 
(Stetson 1992) on the deep, co-added images to obtain instrumental magnitudes of all the 
detected point-sources. The magnitudes of the GRBs were cahbrated relative to bright stars, 
either a single star from the USNO-B catalog or several stars with SDSS photometric data. 
We gathered other R or Rc bands data for these bursts from GCN circulars and the literature 
to complete the optical light curves. While that there are slight differences between R, R^, 
and Sloan r' filters, this has little consequence for our analysis. 

We also reduced the Swift-XKT light curves for the sample. We started from the XRT 
level 2 event files for the photon counting (PC) mode, where events were filtered to be in the 
0.2-10 keV energy band and restricted to grades 0-12. We extracted the XRT-PC spectra 
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and light-curves with the software tool xselect. We used the rmf files from the standard 
XRT calibration distribution, generated the arf files with the Swift-XKT software tool and 
fit the X-ray spectra with XSPEC (Arnaud 1996). 

3. The Afterglow Light Curves 
3.1. GRB 070125 

GRB 0701251 {z = 1.547, Fox et al. 2007) was relatively well observed in the optical 
bands, and GCN data sample between 6 x 10^ and 2.3 x 10^ s after the burst trigger. The 
optical data from 9 x 10^ to 4 x 10^ s can be well fit by a power-law model (dashed line 
in Figure 1) with a slope of a = 1.65 ± 0.05 {x^/dof = 7.4/14). The first LBT data 
point (r' = 26.3 ± 0.3 mag at 26.8 days) deviates by 6.4cr from the single power-law fits. 
The non-detection of the second LBT observation (r' < 26.2 mag at 298 days) indicates 
that the first LBT detection is mostly from the afterglow emission with little host galaxy 
contamination. The data suggests an optical break between the LBT observations and the 
early data taken before 4 x 10^ s, which is also supported by the upper limit of i? > 23.8 mag 
at 10^ s from MDM (Mirabal et al. 2007). We fit the optical data with a broken power-law 
as f{t) oc ((t/t;,)"i* + (t/tb)"^^)'^^^ where we fixed the smoothness parameter to s = 2.5, and 
obtained ai = 1.58±0.12, = 2.87±0.17, and h = (5.0±1.6) x 10^ s with x^/dof = 6.7/13. 

GRB 070125 was also observed by Swift and Chandra. We compared the X-ray light 
curve to the optical light curve, and found that they are broadly consistent with each other. 
We fit the X-ray data before 5 x 10^ s with a power-law to obtain a decay index of ai^x = 1-4- 
1.7 depending on whether we include the X-ray data before 7 x 10'^ s. Although the X- 
ray decay index has a considerable range, even for the extreme case of cti.x = 1-7, the 
extrapolation of the single power-law (4 x 10~^^erg cm~^ s~^) is inconsistent with Chandra 
upper limit of 2 x lO^^^erg cm~^ s^^at 3.4 x 10® s (Cenko et al. 2007), strongly suggesting 
the existence of a break in the X-ray light curve. 

We tested the standard afterglow/jet models (e.g.. Sari et al. 1997; Rhoads 1999 and 
see Zhang & Meszaros 2004 for a review) using our temporal decay indices and the spectral 
index of /3 = 1.03 ± 0.08 that we derived from the X-ray data (consistent with the estimate 
from Burrows & Racusin 2007b). The pre-break spectral and temporal decay indices follow 
the relation ai = 3/?/2 (1.58 ± 0.12 ~ 1.55 ± 0.12), suggesting that the afterglow is at 
the evolutionary stage of < z/ < z/c where and z/c are the typical frequency of the 



^An independent analysis of GRB 070125 is provided by Updike et al. (2008) 
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synchrotron radiation and the coohng frequency, respectively. At this stage, the standard 
models for a constant ISM density predict a post-decay index a2 between 02 = cti + 0.75 ~ 
2.33 ±0.12 and 0:2 = 2/5 + 1 ~ 3.1 ±0.16 depending on whether or not sideways expansion of 
the jet is included. Our post-decay index of 0:2 = 2.87 ± 0.17 is closest to the value expected 
with maximal sideways expansion. 

Thus, GRB 070125 shows a "textbook jet break". Although the optical light curve 
provides much tighter constraints on the pre and post-break decay slopes, the break is 
independently seen in both optical and X-ray light curves, and the consistency of the optical 
and X-ray light curves suggests that the break is achromatic. Furthermore, the spectral and 
temporal decay indices before and after the break satisfy the predictions from the standard 
jet models. 

3.2. GRB 070311 

The late-time optical light curve of GRB 070311 contains a significant fiare at 2 x 10^ s, 
which is immediately followed by a steep decay during the period from 2.3 x 10^ to 3.5 x 10^ s 
with a decay slope of 3.1 ± 0.4 with x^/dof = 0.1/2 (Figure 2 upper panel). A detailed 
analysis of this flare is provided by Guidorzi et al. (2007b). Before the flare, the four optical 
data points (lOMO^ s) set a decay slope of cti = 0.50 ± 0.03 {x^/dof = 4.4/2). After the 
flare, the two LBT data points are significantly below the extrapolation of the power-law 
decay before the fiare, and have a decay slope of 0:2 = 1.6 ± 0.3 {x^/dof = 0/0). This 
indicates a break in the optical light curve, with the break occurring during the fiare from 
10^ s to 3.5 X 10^ s. The post-break decay slope also suggests a relatively low electron index 
of p ~ 1.6. The XRT X-ray light curve also shows a fiare very close to the optical fiare, and 
it follows a similar steep post-fiare decay slope (3.3 ±0.5). The X-ray and optical light curves 
after the fiare and the break are mutually consistent. However, it is difficult to determine 
whether the decay slope for the last few X-ray data points has fiattened to match our LBT 
data given the S/N of the X-ray hght curve. It is also obvious that the optical-to- X-ray 
fiux ratio has changed before and after the fiare/break, suggesting that the early X-ray hght 
curve has a significant fiux contribution from additional emission components. 

3.3. GRB 070411 

GRB 070411 (Figure 2 lower panel) at ^ = 2.954 (Jakobsson et al. 2007) has a well- 
sampled optical hght curve at early times (before 6000 s) and the early light curve shows 
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complicated decay patterns that are also present in many other bursts (e.g., GRB 060206, 
Stanek et al. 2007). The late time optical light curve (after 10^ s), however, is not well 
sampled. We fit the optical light curve after 1000 s with a single power-law and obtained 
a decay index of a = 0.94 ± 0.01 {x^/dof = 119.6/21). Although this power-law fit cannot 
fully represent the complicated early afterglow decay, it is consistent with the late time 
optical light curve, where fitting the data after 4500 s yields a = 0.94 ± 0.02 with improved 
fits with x^/dof — 13.4/6. However, our 3a upper limit from LET at 10.3 days is slightly 
below the extrapolation of the single power-law determined from the early epochs of optical 
data. We also tested fitting the power-law decay to a longer time baseline that included 
the early optical data before 1000 s, and we obtained shghtly flatter slopes which will lead 
to larger deviations from the LET limit at 10.3 days. The X-ray light curve has a slightly 
steeper slope for the single power-law fit {ax = 1.02 ± 0.02 with x^/dof = 19.5/9) than 
the optical slope. Since the last two X-ray data points show possible deviations from the 
single power-law decay, we also tried a broken power-law fit, and obtained a comparable 
fit with x^/dof = 17.9/7). The optical and X-ray data set a lower limit for the jet break 
tb > 2.8 X 10^ s. A break between 3 x 10^ s and 10^ s is consistent with the optical and X-ray 
data; however, the late time data are still not deep enough to exclude the single power-law 
model. 



3.4. GRB 070412, GRB 070419A, and GRB 070518 

GRB 070412 appears optically dark. We observed GRB 070412 twice with the LET at 
2.8 hours and 8.2 days after the burst respectively. We place upper 3a limits of r' > 25.2 
for both observations. Most of the other optical observations of GRB 070412 also only yield 
upper limits. Therefore, we are unable to search for the jet break in GRB 070412. The late 
time optical light curves for both GRB 070419A and GRB 070518 flattened compared to the 
early optical observations taken before the LBT data, indicating other emission components 
besides the afterglow. In GRB 070518, the extra optical flux contribution may come from 
the host galaxy because the LBT images show that the source is extended. GRB 070419A is 
more interesting, since there is a late time bump in the optical light curve, consistent with a 
supernova bump. The details will be discussed in companion papers (Garnavich et al. 2008, 
in preparation; McClelland et al. 2008 in preparation). In both cases, we are unable to 
search for a late time optical jet break. 
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4. Discussion 

We observed six GRBs afterglows using the 8.4m LBT from 2.8 hours to 30.8 days after 
the burst triggers. With 20-30 minute exposures, we are able to reach a flux limit of ~ 25- 
26 mag in the Sloan r' band. These observations systematically probe the late time decay of 
the afterglows. In particular, these observations enable us to extend the temporal baseline 
for measuring the jet breaks by another decade in a logarithmic time scale. 

In three of the six bursts, we are unable to detected the optical jet break because of the 
nature of the afterglow. GRB 070412 is optically dark and the late-time optical light curves 
of GRB 070419A and GRB 070518 contain other emission components from the host galaxy 
or a supernova which will mask a jet break. In the other three cases, we clearly detected jet 
breaks in GRB 070125 and GRB 070311, and a break is probably required for GRB 070411. 
The optical jet breaks in all three cases occur at late-times, when the afterglows have faded 
to be fainter than R = 24 mag. Such faint breaks are not detectable without deep, late-time 
imaging of the afterglow. This strongly suggests that most failures to detect jet breaks are 
due to the lack of well-sampled afterglows either in the rate of sampling or the temporal 
baselines. Several other well-sampled Swift bursts (e.g., GRB 060206, Stanek et al. 2007 and 
GRB 060526, Dai et al. 2007) also clearly show jet breaks. 

The XRT light curves suffer from similar problems. Although the XRT light curves 
have good S/N at early times, the S/N decreases significantly for the time scales on which 
the jet breaks frequently occur. The problem is particularly severe for the last few XRT 
data points, where the count rate is affected by the choice of source region, background 
region, binning, and other effects such as flares either in the GRB itself or in background 
sources. Thus, conclusions about the existence of jet breaks with typical data are heavily 
dependent on the prior for fitting either a single power-law or broken power-law. This is 
also reflected in our sample. For example, in GRB 070125, the XRT data points can be fit 
with both single or broken power-laws. This problem is also present in several other bursts 
(e.g., GRB 060526, Dai et al. 2007), where the optical data require a break while the X-ray 
data are only consistent with a break. In the case of GRB 070125, an upper limit from a 
deep Chandra observation (Cenko et al. 2007) is below the extrapolation from the single 
power-law decay, independently confirming the existence of of a break. 

Analyses of jet breaks are also afi^ected by the complicated decay patterns. When extra 
emission processes, such as energy injection or flares, are important in the light curve, it is 
difficult to determine the forward shock emission decay slope. This may explain the case of 
GRB 070311, where extra emission components in both the optical and X-ray light curves 
may affect the slopes of the afterglows. 
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Given these problems, and the fact that almost all well-sampled bursts in the optical 
bands show a significant break, we argue that the basic picture that GRBs are coUimated 
is still valid. Invalidating the jet model requires examples without jet breaks even at late 
times rather than a failure to find a break at early times. More well-sampled bursts are still 
needed to securely measure a large sample of jet breaks, to constrain the jet models and 
determine the beaming correction for GRBs. Accurate estimates of the beaming corrections 
are needed to determine the intrinsic rate of GRBs and to match those rates to possible 
progenitor populations. 

We would like to thank D. Hartmann, A. Updike, B. Zhang, and N. Gehrels for their 
comments. We thank the Swift team for the prompt detection and localization of GRBs and 
the rapid release of data products, the GRB Coordinates Network (GCN) and astronomers 
who contribute to the GCN circulars. We thank everyone who contributed to the LBT/LBC 
Science Demonstration Time observations. 
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Table 1. LBT Photometry of GRBs 



Burst Observation Days Filter Exposures Mag 

Time After (sec) 

(UT) Trigger 



GRB 070125 


2007-02-21.10 


26.80 


r' 


10 


X 


200 


26.3 ±0.3 




2007-11-19.50 


298.2 


r' 


9 


X 


200 


> 26.2 


GRB 070311 


2007-03-17.14 


6.13 


r' 


10 


X 


200 


24.73 ± 0.06 




2007-03-20.13 


9.05 


r' 


15 


X 


200 


25.42 ±0.12 


GRB 070411 


2007-04-15.15 


3.30 


r' 


10 


X 


200 


23.95 ±0.05 




2007-04-22.14 


10.30 


r' 


10 


X 


200 


> 25.3 


GRB 070412 


2007-04-12.18 


0.12 


r' 


10 


X 


200 


> 25.2 




2007-04-20.31 


8.20 


r' 


5 


X 


200 


> 25.2 


GRB 070419A 


2007-04^23.15 


3.70 


r' 


10 


X 


200 


24.75 ±0.18 




2007-05-10.18 


20.80 


r' 


25 


X 


200 


25.29 ±0.05 




2007-05-20.22 


30.80 


r' 


15 


X 


200 


25.71 ± 0.13 


GRB 070518 


2007-05-19.32 


0.72 


r' 


5 


X 


200 


23.03 ± 0.05 




2007-05-20.31 


1.71 


r' 


5 


X 


200 


23.36 ± 0.05 




2007-05-22.21 


3.60 


r' 


15 


X 


200 


23.56 ± 0.05 
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Fig. 1. — Optical and X-ray light curves of GRB 070125. We included additional R band 
data from the GCN circulars (Cenko & Fox 2007a; Xing et al. 2007; Uemura et al. 2007; 
Haislip et al. 2007; Greco et al. 2007a; Yoshida et al. 2007a; Terra et al. 2007a; Mirabal 
et al. 2007). The dashed line is a single power-law fit (a = 1.65 ± 0.05) to the optical data 
from 9 X 10^ to 4 X 10^ s. The dotted line is a broken power-law fit to the whole optical data 
set with ai = 1.58 ± 0.12, aa = 2.87 ± 0.17, and = (5.0 ± 1.6) x 10^ s. The X-ray and 



optical data are aligned such that the steepest single power fit to the X-ray data {ax 
overlaps the single power fit to the optical data (the dashed line). 
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Fig. 2.— Optical and X-ray light curves of GRB 070311 and GRB 070411. We included 
additional R band data from Guidorzi et al. (2007) and GCN circulars (Yosliida et al. 
2007b; Halpern & Armstrong 2007abcd; Jelinek et al. 2007a; Updike et al. 2007a; Greco et 
al. 2007b; Kann et al. 2007a) for GRB 070311, and for GRB 070411 (Mikuz et al. 2007a; 
Gomboc et al. 2007; Jelinek et al. 2007b; Berger et al. 2007; Mikuz et al. 2007b; Kann et 
al. 2007b; Ferrero et al. 2007; Perley et al. 2007). 



